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Abstract 

In this paper, by both simulations and theoretical predictions we study two and 
three node (or degree) correlations in random Apollonian network (RAN), which 
have small-world and scale-free topologies. Using the rate equation approach under 
the assumption of continuous degree, we first give the analytical solution for two 
node correlations, expressed by average nearest-neighbor degree (ANND). Then, we 
revisit the degree distribution of RAN using rate equation method and get the exact 
connection distribution, based on which we derive a more accurate result for mean 
clustering coefficient as an average quantity of three degree correlations than the 
one previously reported. Analytical results reveal that ANND has no correlations 
with respect to degree, while clustering coefficient is dependent on degree, showing 
a power-law behavior as C{k) ~ k^ 1 . The obtained expressions are successfully 
contrasted with extensive numerical simulations. 

Key words: Complex networks, Correlations, Scale-free networks, Small-world 
networks, Networks 



1 Introduction 



The last few years have witnessed the birth of a new movement of interest 
and research in the study of complex networks as an interdisciplinary subject 
[1,2,3,4,5]. This flurry of activity, triggered by two ground-breaking papers by 
Watts and Strogatz on small-world networks [6] and Barabasi and Albert on 
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scale- free networks [7], has been certainly induced by the increased comput- 
ing powers and by the possibility to study the properties of a plenty of large 
databases of real-life networks, such as Internet, World Wide Web, metabolic 
networks, protein networks in the cell, co-author networks and citation net- 
works. These networks have been identified and analyzed in the literature, 
the emphasis being mostly on the three basic topological characteristics, such 
as power-law degree distribution, small average path length (APL) and high 
clustering coefficient [1,2,3,4,5]. 

However, the above mentioned three properties do not provide sufficient char- 
acterizations of the real- world systems. In fact, it has been observed that real 
networks exhibit ubiquitous degree correlations among their nodes [8,9,10,11,12] 
This translates in the observation that the degrees of nearest neighbor nodes 
are not statistically independent but mutually correlated in most real-life net- 
works. Correlations play an important role in the characterization of net- 
work topology, and have led to a first classification of complex networks [11]. 
They are thus a very relevant issue, especially in view of the important con- 
sequences that they can have on dynamical processes taking place on net- 
works [13,14,15,16]. 

Recently, based on the well-known Apollonian packing, deterministic Apollo- 
nian network (DAN) [17] and random Apollonian network (RAN) [18], called 
jointly Apollonian networks (ANs) were proposed, which have been general- 
ized to high dimension [19,20,21,22]. ANs are simultaneously scale-free, small- 
world, Euclidean (i.e. it can be embedded in an Euclidean lattice [17]), space 
filling, and with matching graphs (that is to say, in the infinite network size 
limit, the edges do not only completely cover the space like a Peano curve, 
but also never cross each other). They may provide valuable insight into the 
real-life networks, e.g. the maximal planarity [18] of ANs is of particular prac- 
ticability for the layout of printed circuits. Moreover, ANs describe force chains 
in polydisperse granular packings and could also be applied to the geometry of 
fully fragmented porous media, hierarchical road systems, and area-covering 
electrical supply networks [17]. Particularly, DAN can be helpful to understand 
the energy landscape networks [19,23,24]. Very recently, ANs have attracted 
increasing interest from the scientific community [25,26,27,28,29,30], some in- 
teresting dynamical processes, such as percolation, epidemic spreading, syn- 
chronization and random walks taking place on RAN have been investigated. 

In this paper we investigate the correlation properties of nodes' connectivity of 
random Apollonian network (RAN) [18,21,22]. Combining a rate equation [31] 
in the continuous degree approximation and the boundary condition [32] of 
this rate equation, we work out analytically the two node correlations in RAN, 
measured by the average nearest-neighbor degree (ANND) of nodes with de- 
gree k. Then, using the genuine discrete degree distribution, we obtain the 
exact analytical result of mean clustering coefficient, which is the average 
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value of three node correlations. Both the analytical expressions are in very 
good agreement with numerical simulations. 



2 Random Apollonian network 



The random Apollonian network (RAN) [18] is the stochastic version of the 
deterministic Apollonian network (DAN) [17]. The construction of DAN begins 
with (t — 0) a triangle and at each time step a new node is connected to the 
nodes of every existing triangles, omitting those triangles that had already 
been updated in previous steps. In the random version the triangles to be 
updated are selected at random, one at a time. We call the triangles to be 
updated active triangles. 

Here we focus on the random Apollonian network (RAN). Some properties 
of the RAN have been investigated [18,21,22]. The average degree of all its 
nodes equals 6. It has a power-law degree distribution P(k) ~ k~ 3 , and since 
each new node induces the addition of three triangles, one expects that the 
RAN has a finite clustering coefficient. It has been proved that the increasing 
tendency of average path length of RAN is a little slower than the logarithm 
of the number of nodes in RAN. 

For the sake of the following investigation on correlations in RAN, we first 
revisit the degree distribution of RAN. Since the network size is incremented 
by one with each time step, we use the step value t to represent a node created 
at this step. Note that at time t there are 2t + 1 active triangles in RAN, and 
the number of active triangles containing a node with degree k is also equal 
to k. Let N k {t) denote the average number of nodes with degree k at time t. 
When a new node enters the network, N k (t) changes as [33] 

alN k (k - l)N k _i(t) - kN k (t) 

= r k 1 

alt 2t+l fe ' 3 K J 

In the asymptotic limit N k (t) = tP(k), where P(k) is the degree distribution. 
Eq. (1) leads to the following recursive equation 



P(k) = < 



gP(fc - 1) for k > 3 + 1 



( 2 ) 



for k = 3 , 



giving 

24 

P{k) = k(k + l)(k + 2Y (3) 

In the limit of large k, P(k) ~ A; -3 [18,21], which has the same degree exponent 
as the BA model [7] and the hierarchical lattice [34,35]. This obtained degree 
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distribution is a more accurate result than the previous one [18,21]. 



Notice that the growing precess of RAN actually contains the preferential 
attachment mechanism, which arises in it not because of some special rule in- 
cluding a function of degree as in Ref. [7] but naturally. Indeed, the probability 
that a new node created at time t will be connected to an existing node % is 
clearly proportional to the number of active triangles containing i, i.e. to its 
degree ki(t). Thus a node % is selected with the usual preferential attachment 
probability IL[/cj(t)] = ki(t)/(2t+ 1) ~ ki(t)/2t (for large t). Consequently, /c, 
satisfies the dynamical equation [1] 

dt 2t ' 1 ' 

Considering the initial condition ki(i) = 3, we have [21] 

m 1/2 

Ht) = 3 y . (5) 

Equation (5) shows that the degree of all nodes evolves the same way, following 
a power law as in the well-known BA networks [7]. 



3 Correlations in Random Apollonian network 



Having obtained the exact degree distribution of the random Apollonian net- 
works (RAN), we now study the two and three node correlations in RAN, 
which are merely two specific types of correlations. We ignore the long-range 
and multinode correlations, since the empirical data and theoretic research on 
such correlations is also absent. 



3.1 Two node correlations 



Two node correlations in a network can be conveniently measured by means 
of the quantity, called average nearest-neighbor degree (ANND), which is a 
function of node degree, and is more convenient and practical in characterizing 
degree-degree correlations. The ANND is defined by [9] 

k nn (k) = Y / k'P(k'\k). (6) 

k' 

If there are no two degree correlations, k nn (k) is independent of k. When 
k nn (k) increases (or decreases) with k, the network is is said to be assortative 
(or disassortative) [11,12]. 
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We can analytically compute the function value of k nn (k) for the RAN. Let 
Ri(t) denote the sum of the degrees of the neighbors of node i, evaluated at 
time t. It is represented as 

R t (t) = £ kj(t), (7) 
jev(i) 

where V(i) corresponds to the set of neighbors of node i. The average degree 
of nearest neighbors of node i at time t, k nn (i,t), is then given by k nn {i,t) = 
Ri{t)/ki(t). During the growth of the RAN, Ri(t) can only increase by the 
addition of a new node connected either directly to i, or to one of the neighbors 
of i. In the first case Ri(t) increases by 3 (the degree of the newly created 
node), while in the second case it increases by one unit. Therefore, in the 
continuous k approximation [1,2,3,4,5], we can write down the following rate 
equation [31,32]: 



d -^ = m l [k l {t)]+ E n#,(t)] = ^ + ^. (8) 

ie(V)(i) 

The general solution of Eq. (8) is 

i 9 n\ 1/2 

R i (t) = Qo(i)t* + -[-) \nt, (9) 



2 \i 

where $o(i) is determined by the boundary condition Ri{i). To obtain the 
boundary condition Ri(i), we observe that at time i, the new node i is con- 
nected to an existing node j of degree kj{i) with probability Uj[kj(i)}, and 
that the degree of this node increase by one unit in the process. Thus, 



R l (t) = T, Il AkAm^) + M- (io) 

1/2 

Inserting Tlj[kj(i)} = and kj(i) = 3 (jj into Ri(i) leads to 



^« = 3 + -E-^3 + -ln^. (11) 
j=i J 

So, in the large i limit, Ri(i) is dominated by the second term, yielding 

Ri(i) = ^hii. (12) 

From here, we have 

Ri(t)*l(j) ht, (13) 
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Fig. 1. Plot of average nearest-neighbor degree of the nodes with degree k. The 
squares denote the simulation results for various network sizes, while the solid lines 
are the theoretical result provided by Eq. (14). 



and finally 

k nn (k,t) ~ -Int. (14) 

So, two node correlations do not depend on the degree. The ANND grows 
with the network size N = t as hat, in the same way as in the BA model [10]. 

In order to confirm the validity of the obtained analytical prediction of ANND, 
we performed extensive numerical simulations of the RAN (see Fig. 1). To 
reduce the effect of fluctuation on simulation results, the simulation results 
are average over fifty network realizations. From Fig. 1 we observe that for 
large k the ANND of numerical and analytical results are in agreement with 
each other, while the simulated results of ANND of small k have a very weak 
dependence on k, which is similar to the phenomena observed in the BA 
model [10]. This k dependence, for small degree, cannot be detected by rate 
equation approach, since it has been formulated in the continuous degree k 
approximation. 

From Eq. (14), we can also easily know that the evolution of ANND of all 
nodes is the same, showing a weak logarithmic divergence with evolution time 
(network size). To confirm this analytical prediction we also perform numerical 
simulations. In Fig. 2 we report the numerical results, which are in good 
agreement with this prediction provided by Eq. (14). 
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Fig. 2. Semilogarithmic graph of time evolution for the ANND of the node added 
to the system at time 2. The straight line follows k nn (k,t) = | hit as predicted by 
Eq. (14). All data are from the average of 50 independent runs. 

3.2 Three node correlations 



Three node correlations can be measured by means of the conditional proba- 
bility P(k', k"\k) that a node of degree k is simultaneously connected to nodes 
with degree k' and k" . The direct evaluation of the P(k',k"\k) is generally 
difficult. To overcome this problem, another interesting alternative quantity, 
i.e. clustering coefficient, is frequently used. By definition, the clustering coef- 
ficient [6] of nodes with degree k, C(k), is defined as the probability that two 
neighbors of a node with degree k are also neighbors themselves. The average 
clustering coefficient of the whole network is then given as 

C = Y,P(k)C(k). (15) 

k 

Now we analytically estimate the average clustering coefficient C of RAN by 
means of the clustering spectrum C(k). In the RAN, for a node of degree k, 
the exact value of its clustering coefficient is [18,21,22] 

4k — 6 

c < fc > = W^ry < 16 > 

which depends on degree k. For large k, the degree-dependent clustering C(k) 
is inversely proportional to degree k, the same behavior has been analytically 
found in some growing network models such as deterministic (random) pseud- 
ofractal scale-free networks and their variants [36,37,38,39,40,41,42], as well 
as some real systems [43]. 

Using Eq. (15), the average clustering coefficient C of RAN can be easily 
obtained as the mean value of C{k) with respect to the degree distribution 
P{k) expressed by Eq. (3). The result is 
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k=3 k(k + l)(k + 2)k(k-l) 

fz^-WA-^W^, (it) 



fc=3 



+ l)(fc + 2) \k-l k(k-l) 



where d = Eg* fc(fc+1 2 } 4 (fc+2) ^ and C 2 = ESf =3 ^WAv We now 
compute in detail Ci and C 2 , respectively. First, we can decompose C\ into 
the sum of four terms as 



~ / 16 48 48 16 \ 4 

Cl " ^ 3 life^i ~T + a 7 TT~ fcT2 J " 3" (18) 

Analogously to Eq. (18), we get 



~ / 24 36 72 72 12 \ 9 



fc=3 



where we have used the fact that Em=i ^2 = g^2- Substituting Eqs. (18) and 
(19) into Eq. (17), we have 



4 / \ 
C =- - (119 - 127T 2 ) ~ 0.768 

Thus the average clustering coefficient C of RAN is large and independent of 
network size. Since Eqs. (3), (16), (18), and (19) are exact, it is the same with 
the obtained C value. We have performed extensive numerical simulations of 
the RAN. In Fig. 3, we present the simulation results about the average clus- 
tering coefficient of RAN, which are in complete agreement with the analytical 
value. 



(20) 



It should be mentioned that previous jobs [18,21] have studied the average 
clustering coefficient of RAN, using continuum approximation and integral 
methods, and have acquired a analytical value of C = y — 361n| ~ 0.737. 
From Fig. 3, we can see that there is a difference between this analytical 
value and simulated ones. Where does the difference come from? And why Eq. 
(20) is more accurate? These might be explained as follows. In RAN, a node 
may have only an integer number of connections (degrees). The continuum 
approach used previously [18,21] is under the assumption that the degree 
distribution P(k) is modeled by the continuous distribution P c (k) = ak~ 3 (k > 
3), where a = 18 is a normalization constant. It does not properly account 
for the fraction (i.e. degree distribution) of nodes with small degrees [44], and 
therefore has a discrepancy with the exact result for RAN with the genuine 
discrete degree distribution Pd(k) expressed by Eq. (3) [see Fig. 3]. However, 
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Fig. 3. Average clustering coefficient C of RAN vs the network size N. The solid line 
shows the exact result | — (119 — 12-7T 2 ), while the dashed line depicts the previous 
analytic value — 36 In |. The squares denote the results of simulations. Each data 
is obtained by fifty independent network realizations. 

it is worth noticing that (a) P c {k) and Pd{k) are equal, asymptotically, in the 
limit of large k and (b) the difference is most pronounced for k w 3. The 
discrepancy is mainly due to the values of P(k) for small k. Thus, the discrete 
method, omitting a continuum assumption, seems more suitable for obtaining 
exact results for some challenging tasks [1]. 



4 Conclusion 

In this article, we have done an analytical study of the correlations between 
the degrees of neighboring nodes in random Apollonian networks (RAN). Ap- 
plying rate equation method in the continuous k approximation, together with 
its boundary condition that node % was added to the system at time i with 
the expected sum of degrees of its neighbors being Ri(i) =3 + 1 lnz, we have 
provided the solution of the average degree of nearest neighbors as the mea- 
sure of two degree correlations. The obtained result shows that RAN lacks 
two degree correlations, which means that the ANND of nodes with degree 
k is independent of k. Moreover, we have obtained the exact degree distri- 
bution of RAN, on the basis of which we have attained the exact value of 
mean clustering coefficient as an average quantity of three node correlations. 
We found that RAN are highly clustered, and the clustering spectrum scales 
as C(k) ~ k~ l , exhibiting hierarchical organization and modularity. These 
obtained structure properties may be helpful to understand and explain the 
workings of systems built on RAN. Many physical processes, such as disease 
spread and random walks, taking place in RAN exhibit different behaviors 
than those in the classic BA networks, which are relevant to these particular 
topologies of RAN [18,26,30]. 
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